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1. Introduction 

Lattice calculations allow us to study hadronic properties from first principles. One of the 
great challenges of lattice QCD, however, is that the physical light quark masses, in particular the 
physical up and down sea quark masses, are numerically inaccessible with present hardware. We 
thus rely on (partially quenched) chiral extrapolations to get to the physical light quark regime 
from simulations with unphysically heavy quarks. We can also determine the low energy constants 
(LECs) [|l|], aka Gasser-Leutwyler coefficients, which appear in the chiral Lagrangian, from these 
partially quenched simulations provided the systematic errors are treated properly and the correct 
number of fermion flavors are used [^]. 

Domain wall fermions (DWF) have the advantages of nearly exact chiral symmetries with a 
moderate extent in the fifth dimension and a continuum-like chiral perturbation theory. The resid- 
ual chiral symmetry breaking of DWF, quantified as mres [||], comes in as an additive constant to 
the quark masses at the lowest order approximation. The employment of the rational hybrid Monte 
Carlo (RHMC) technique [Q], free of finite step-size errors, also gives us better control over system- 
atic errors. It thus makes possible a direct comparison of our numerical results with the predictions 
of chiral perturbation theory. This proceeding focuses on a set of A/j = 2 + 1 DWF simulations on 
16'' X 32 X 8 lattices with a lattice spacing of 1.6(1) GeV. After detailing the simulation parameters, 
I present results for the masses and decay constants of the pseudoscalar mesons. Preliminary chiral 
fits up to next-to-leading order (NLO) on these quantities are then given. A rough investigation of 
the effect of nondegeneracy is also described. 

2. Simulation details 

The gauge action used in our simulation can be written in the general form : 

5G[t/] = -f ((l-8ci) £ P[[/],,^v + ci £ /?[f/k^.v) (2.1) 

We generated three sets of 16^ x 32 x 8 gauge configurations with the DBW2 action, in which 
c\ = —1.4069, at j3 = 0.72. Two light sea quarks with equal mass m/ and one strange quark with 
mass mv were included in the fermion determinant. The RHMC algorithm was used to update 
the gauge fields. All the data was generated on QCDOC machines at Columbia, the RBRC and 
Edinburgh. 

Valence measurements with two degenerate quarks of up to 8 different masses my were per- 
formed on these lattices. Nondegenerate valence measurements with light quark masses of 0.005, 
0.01, 0.015, 0.02 and a heavy quark mass of 0.04 were also done for mi = 0.01 and 0.02. We 
skipped the first 1000 trajectories for thermalization and measured point-point meson correlators 
thereafter. Table [T] shows the relevant parameters for evolution and measurements. The integrated 
autocorrelation time was determined to be on the order of 50 to 100 trajectories. Thus we binned 
the data into blocks of 100 trajectories prior to statistical analysis to account for the autocorrelations 
in a robust fashion. All the quantities are in lattice units unless noted. 

The residual chiral symmetry breaking turns out to be quite large at this coupling. When 
extrapolated to the dynamical chiral limit (nii = 0), 

mres = 0.0106(1) (2.2) 
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Table 1: Simulation parameters. 5t is the step size in molecular dynamics integration and T is the total 
length in one trajectory. Data for the m/ = 0.01 and 0.02 ensembles used in this proceeding reflects twice as 
much statistics as that in the talk given at the conference. 
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which is comparable to the input light sea quark mass. Noting that the coupling is fairly large and 
Ls is rather small, it is not surprising to have a large wires- For hadronic observables like meson 
masses and decay constants, we treat mres as a shift to the input quark masses and neglect other 
possible higher-order effects. 



3. Preliminary Results 

The pseudoscalar masses and decay constants can be extracted from both the pseudoscalar 
and axialvector correlators [||]. Here emphasis will be given to the results from the pseudoscalar 
channel, since they give smaller statistical errors than the axialvector correlators. 

Mji and the chiral fits. The next-to-leading order (NLO) quark mass dependence of the pseu- 
doscalar masses in PQ;^PT with N flavors of sea quark has been computed generally [^]. With 2-1-1 
flavors of sea quark and two degenerate valence quarks, the formula simplifies to 

Ml = xv[y+^^ (2L6 - + ^ (2^8 - Ls)Xv 

1 \ ^Xv-Xi-Xs ^, {Xv - Xi){Xv - Xs) ^, 

+ - ^'X^" - (■ + iog;tv)+ ^-' z,iog;t,] } (3.1) 

Xv Xr\ yXv Xr\) ^ 

where Xx = 2Bm^, x = V, I, s, or t], Ztj = + ^Z^) and x = ^{'^Xi + Xs)- Here we have 
explicitly taken the chiral scale to be 1 GeV. For DWF, the masses should all be shifted by ntres- 
Thus to leading order in the valence quark masses, = 2B{mv + mres)- 
As a consistency check, we fit our results for to the linear form 

Ml = 2B{mv + mres) + C (3.2) 

where mres is the residual mass at the valence chiral limit (my — > 0), i.e., 0.0111(1), 0.0113(1) and 
0.0122(1) for m/ = 0.01, 0.02 and 0.04 respectively. We would expect C to vanish if the lowest- 
order approximation is good enough. Independent linear fits for the pion masses with different sea 
quarks are shown in the left panel of Figure [l|. The two heaviest masses were excluded from the 
fit except for m/ = 0.04 where all the available data were used. In all the cases, is quite close 
to zero when my = —m^^^. The larger deviation from for m/ = 0.02 may be due to higher order 
corrections from ;^PT or ^(a^ ) violations of chiral symmetry, which need to be further investigated. 
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Figure 1: LEFT: vs. my in the pseudoscalar channel; RIGHT: M'^/{mv +»Jres) vs. my. The lines 
are fits to the mass range [0.005,0.03]. Filled symbols and open symbols represent data from measurements 
with degenerate and nondegenerate quarks, respectively, my is the average of two valence quark masses. 



We also show M|/ (my + mres) as a function of niy in the right panel of Figure Nonlinearities 
are evident in all the datasets. The solid lines are the partially quenched NLO fit to Eq. 3.1 with 
sea quark mass fixed. We chose the fitting range to be [0.005,0.03]. In the fit, the value of at 
mi = my = — '«res from the dynamical linear extrapolation was used as the input for /, while in 
principle the limit of lUs —wires should also be taken to obtain the value for /. As can be seen, 
for mi = 0.01 the curve does go through all the data. For mi = 0.02 the NLO fit represents the data 
well, for all but the heaviest valence masses, where differences of about one standard deviation are 
seen. As will be shown later, the NLO fit to fj^ also has the same problem. 

fjt and the chiral fits. We use the definition of fji in which its physical value is about 130 
IVIeV. When two degenerate valence quarks are used, the NLO mass dependence of fji can be 
simplified to ^ 

h = f[ 1 + J2{NLaX+L5Xv) - ^^^2f2 [iXv+Xi)log — - — + — - — log — - — J j(3.3) 

We show the results for /;;: from the pseudoscalar correlators in Figure ^. The data points have 
obvious curvatures, especially for mi = 0.02. The straight lines fit the data poorly. We also fit the 
mi = 0.01 and mi = 0.02 data independently to Eq. 33 using the values of B from the NLO fit 
of Ml as an input. Two fitting ranges were used: [0.005,0.025] and [0.005,0.03]. For m/ = 0.01, 
the fitting results do not depend on the fitting range we choose. Even a fit to the first five points 
reproduces the data at heavier quark masses quite well. However, this is not the case for m/ = 0.02, 
where the fit curve misses the data points at heavier masses. It is especially obvious for the fit to 
the range [0.005,0.025]. Similar behaviors were observed in Nf = 2 DWF simulations [^], where, 
instead of fitting data with different sea quarks independently, simultaneous fits to fj^ from different 
sea quarks were performed. 

Remarks. We cannot conclude here that we have found a signal of chiral logs. The nonlin- 
earities we see for both M| and /;;: may also be attributed to finite-volume corrections, zero-mode 
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Figure 2: f„ with linear and NLO fits to mass range [0.005,0.025] and [0.005,0.03] for m, = 0.01 (left) 
and mi = 0.02 (right) datasets. Filled symbols and open symbols represent data from measurements with 
degenerate and nondegenerate quarks, respectively, my is the average of two valence quark masses. 

effects or uncertainties in the residual mass. For the large quark masses we used, NNLO contribu- 
tions may already be important, especially for /^j. More thorough investigations of these systematic 
effects are needed to justify our NLO fits. 

4. Effect of nondegeneracy 

We also explored the possible effects of nondegenerate quarks in the pseudoscalar masses and 
decay constants. Suppose we have two sets of valence measurements on the same lattices. One uses 
two degenerate quarks (m^ = mg = my) and the other uses two nondegenerate quarks (m^ 7^ ms, 
niy = Leading order ;^PT predicts, for the pseudoscalar masses, = B{mA + mB). If we 

take the nondegenerate dynamical points (m^ = mi, = rris) and the matching degenerate points 
(my + mres = ) to extrapolate to the light sea quark limit (m; —>■ — mres), we would expect 

the results to be the same assuming the higher order corrections are small. This gives us rough 
estimates of the kaon mass from both the degenerate and nondegenerate valence quarks. The same 
comparison can also be done for the pseudoscalar decay constants. We show the comparisons in 
Figure ||. The two extrapolations gave quite consistent results within errors. 

5. Summary and Outlook 

Our present data are still statistically limited. Further improvements are needed to reduce the 
statistical uncertainties in the determinations of Mj^ and /;;:. In the meanwhile, better estimates of 
the systematic effects coming from finite volume, finite lattice spacing and residual chiral symmetry 
breaking are necessary to justify our chiral fits and determine LECs. For the work reported here, 
the smallest Mj^L is 4.2, which should make the finite volume corrections few percent effects. The 
ongoing 24^ x 64 x 16 simulations with 2+1 dynamical flavors of DWF on QCDOC will provide 
us better opportunities to probe the chiral limit of various physical quantities. 
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Figure 3: Estimate the effects of nondegeneracy. 
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